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Abstract 

Self-similar solutions of the coherent diffusion equation are derived and measured. The set of 
real similarity solutions is generalized by the introduction of a nonuniform phase surface, based 
on the elegant Gaussian modes of optical diffraction. In an experiment of light storage in a gas of 
diffusing atoms, a complex initial condition is imprinted, and its diffusion dynamics is monitored. 
The self-similarity of both the amplitude and the phase pattern is demonstrated, and an algebraic 
decay associated with the mode order is measured. Notably, as opposed to a regular diffusion 
spreading, a self-similar contraction of a special subset of the solutions is predicted and observed. 
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Self-similar solutions are generally associated with the long-time behavior of dynamic 
processes [l] and found in nearly all disciplines, from astrophysics and fluid dynamics to 
condensed matter and optics . In most dissipative systems, similarity solutions decay 
with a characteristic rate, indicating the asymptotic evolution of a given initial condition 
6|, |8| . Similarity solutions emerge also in non-dissipative systems and often prevail, with a 
familiar example being the family of Gaussian beams in free-space paraxial optics [9| , when 
the propagation distance is given the role of time. Gaussian beams are broadly referred to 
as modes of diffraction, even though they are not nondiffracting as Bessel beams 10] and 
not eigenmodes of an underlying Hamiltonian. Moreover, they are only partially self-similar 
throughout the propagation - their shape is preserved up to scaling and normalization, while 
their phase pattern curves or flattens. 

The imaginary-time counterpart of paraxial diffraction is the coherent diffusion of a 
complex-valued field, 

^ = DVli/i - 70, (1) 

with ip = ifj(x,y;t) being a two-dimensional coherence field, D a real coefficient, and 7 
a linear decay rate. Coherent diffusion in the form of Eq. (DO) arise, for example, in the 
thermal motion of hot coherent atoms, where the field of internal-state quantum coherence 

n 

(ip) is subjected to atomic diffusion (D) and decoherence (7) [111 ]. Using the technique 
of light storage and retrieval [ljj], any arbitrary initial condition can be imprinted on the 



diffusing atoms, and the subsequent dynamics can be observed 



131 ] . This system attracted 



considerable recent study, exhibiting spectral fringes, narrowing, and coherent recurrence 
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1611 . magnetization diffusion 17f, diffusion of vortices la l. transverse momentum diffu- 

fin n 

sion [19|, |20|], and slow-light spreading [211 ] . 

In this letter, we present and explore the exact self-similar complex solutions of Eq. ([Q). 
We experimentally create and follow the dynamics of several self-similar modes and demon- 
strate the preservation of both their shape and phase pattern, as well as their characteristic 
decay. A self-similar contraction, resembling focus or collapse ]2], is also demonstrated. 

Consider first the time- independent paraxial diffraction, dE/dz = —iV 2 j_E/(2k), for the 
slowly-varying envelope, E, of a light field with a wave number k. Two different sets of 
polynomial-Gaussian solutions are known for this equation, namely the standard and the 
elegant beams (9]. The more familiar 'standard' modes, e.g. the Hermite-Gaussian (HG) or 
Laguerre-Gaussian (LG), form complete sets of modes that are self-similar under diffraction. 




FIG. 1: Top: the effect of diffusion on a Gaussian beam is an effective stretching of the waist 
radius, even when the diffusion does not occur at the waist plane. Therefore, far enough from 
the waist, the diffusion results in the contraction of the transverse shape. Bottom: colormap of 
the phase gradients of the field envelope (black lines are equal phase contours; white lines are the 
beam outline). From the viewpoint of the microscopic atomic motion, the contraction far from 
the waist occurs due to destructive interference of atoms diffusing through the rapidly oscillating 
phase pattern (red colored). 



Their transverse intensity distribution, / (x, y; z), is maintained along the propagation direc- 
tion z, normalized, and scaled by the beam radius w(z). In contrast, the transverse shape 
of the 'elegant' solutions is generally not maintained, and originally they were investigated 
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due to their elegant mathematical form 

The elegant HG solution, E^ n2 (r; w ), with w being the radius at the waist plane z = 0, 



is written in terms of the Hermite polynomials of orders n\ and ri2 as 
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Here, q(z;wo) = %zr + z is the complex radius, zr = kw\j1 is the Rayleigh length, iV = 
n i + n 2 is the total mode order, and E is a normalization constant. The transverse scaling, 
appearing in both the polynomial and the Gaussian terms, depends on z and wq, with 
x = x{ik/2/q(z; Wo)] 1 / 2 and y = y{ik/2/q(z;w )} 1 / 2 . The beam radius, w(z), and the radius 
of curvature of the phase fronts, R(z), are obtained from q(z; Wq)^ 1 = R(z)^ 1 — (i2/k)w(z)~ 2 . 
For the corresponding standard mode, the complex arguments of the polynomial (x, y) are 
replaced by real arguments [y/2x/w(z), \/2y / w{z)]. Thus elegant modes with a homogeneous 
polynomial are also standard, and we denote them as 'common' modes. 
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Now suppose that a two-dimensional diffusion takes place at a certain (x, y) plane, where 
z is held constant. Eq. (J2J), with dE/dz = —i\7 2 ± E/(2k) and (d/dw ) z = ikw {d/dq) Wo + 
(d/dw ) g , then gives 

2 / dE RG \ 2 
Vi<? Ba = - " A ( 3 ) 

Under diffusion, the first term in the right-hand side of Eq. ([3]) accounts for a change in the 
waist radius, and the second term for an homogenous decay of the field. Eq. (CQ) is therefore 
solved by 

(x, y, t) = e-^s (t)- (JV+1) E™ n2 (r; w s (t)) , (4) 

where the diffusion coefficient enters only through the waist stretching factor, s (t) = (1 + 
ADt/wl) 1/2 . Thus the spatial consequence of diffusion is always an effective stretching of the 
beam radius at the waist plane, even when the diffusion occurs far from the waist (z ^ 0), 
as illustrated in Fig. [T] (top). Note that if diffusion is addressed as an imaginary-time 
evolution of diffraction, it is readily seen from the definition of the complex radius q(z; wq) 
that exchanging the real evolution in z for an imaginary one corresponds to a real increase in 
the waist radius. The total power, P(t) = J dxdy\ip^^ z 2 \ 2 , which is independent of z, is not 
preserved under diffusion even when 7 = 0, due to an algebraic decay term s (t)~^ N+1 \ This 
occurs even for the lowest Gaussian mode (N = 0) because it is the field, rather than the 
intensity, that is diffusing. Higher-order modes (N > 0) decay faster due to the diffusion of 
the non-homogenous phase pattern, which contains larger gradients for higher N. A similar 
procedure can be carried out for an elegant LG solution, E^, of radial order p and orbital 
order m, yielding Eq. (TJJ for the diffusing field ^ ,z (s, y, t), with N = p + m. 

At the waist, all arguments in Eq. ([2]) are real, and the HG solutions ^ =0 are iden- 
tified with the expanding similarity solutions of regular diffusion [8], occurring, e.g., for the 
vorticity field of a viscous fluid 6|, |24j]. The real solutions are alternatively derived from 
a given self-similar solution with a single real scaling (ip(r,t) = h(t)f[r/w(t)]), by taking 
any of its spatial derivatives (d n ip/dx n = h(t)f^[r/w(t)]/w(t) n ). Indeed, the derivatives 
of the lowest-order Gaussian beam constitute the elegant modes, with the derivative order 
corresponding to the total mode order N, and with the possible generalization for unified 
Hermite-Laguerre-Gaussian modes [25|, [26|. Consequently, ipn^n 2 = ° or V ; p,m Z= °; an d a U linear 
combinations of them with the same total order, are self-similar, sharing the same stretching 



4 





FIG. 2: (a) Diffusion of the complex field of atomic coherence in a storage-of-light experiment, 
demonstrating the self-similarity property of (top to bottom): the basic Gaussian mode, the 
Laguerre-Gaussian modes with radial order p = and orbital orders m = 1,2, and the Hermite- 
Gaussian mode with cartesian orders (0, 1). All images are 1.6 x 1.6 mm. (b) Images after storage, 
diffracted by a binary grating mask with a fork dislocation, for confirming the conservation of phase. 
The retrieved Gaussian mode (top) yields the two m = ±1 vortex modes, while the retrieved vortex 
m = +1 (bottom) produces an m = and an m = +2 modes. 



and the same algebraic decay. Alternatively, any complex 'image' in the (x, y) plane can be 
expanded in terms of ipnun 2 = ° or ^p,™^ using their biorthogonal pairs [27j , its diffusion can 
be described in terms of the modes dynamics, and asymptotically the lowest-order solution 
prevails [30]. 

To validate the above predictions, we perform experiments with thermal alkali atoms 
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confined in a vapor cell with a buffer gas. A resonant laser beam with the desired optical 
mode is sent into the cell, and its complex field envelope is mapped onto the atomic coherence 
field utilizing a storage-of- light technique, by shutting down an auxiliary control beam 
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131 ] . The coherence field, ip, is allowed to evolve for a controlled duration r, in which the 
alkali diffusion through the buffer gas takes place, as well as a homogenous decoherence (e.g. 
due to spin-exchange relaxation 28|). The coherence is then converted back to light, which 
is imaged onto a camera. The experimental set-up and procedure is similar to that described 
where the topological stability of the vortex in a stored Eq± mode was attested. 



m 



The experiment was carried out with the fundamental Gaussian mode E^^^ LG \ the LG 
modes E^f and Eq%, and the HG mode E^f . Fig. |2(a) presents the retrieved images, 
proportional to \ip(x, y, r)| 2 , for a storage performed with the cell located at the beam waist 
(z — 0), for durations of r = 2, 30, and 60 [is. Evidently, all the modes expand but maintain 
their shape through the diffusion process. As a complementary test, we have also passed 
the retrieved beams through a binary grating mask with a fork dislocation, which adds a 
phase function mcj) in its m-th diffraction order 29]. After the mask, as shown in Fig. EJ^b), 
the retrieved vortex mode E^f (m = 1) produces a Gaussian (m = 0) and a higher-order 
vortex (m = 2) in the —1 and +1 diffraction orders, confirming the maintenance of the 
phase pattern. Fig. |3]^a) presents the increase in the waist-radii squared versus the storage 
duration, showing the same linear increase for all curves, w(r) 2 — w$ = Wq[s (t) 2 — 1] = 4Dr. 
The cross-sections shown in Fig. [3](b) , scaled according to s (r) and normalized, clearly 
demonstrate the self-similarity. The algebraic decay of the diffusing modes, s (t) ^ 4 ^ 1 ^ , is 
measured by integrating over the intensity of the retrieved images (Fig. Hj). All modes 
exhibit a significant algebraic decay on top of the homogenous decay, with the higher-order 
modes decaying faster, showing an excellent quantitative agreement with the predictions of 
Eq.©. 

We now discuss the diffusion of the elegant modes at a given plane outside the waist 
plane, z ^ 0. For Gaussian modes, an expansion of the waist radius results in an increase 
of the beam's transverse size for \z\ < zr and in a decrease for |z| > zr (Fig. (TJ top). It 
follows, perhaps counterintuitively, that the initial effect of diffusion occurring at \z\ > zr 
is a contraction. Locally, it is the consequence of a destructive interference in regions of the 
beam where the phase pattern rapidly oscillates (Fig. [TJ bottom). The contraction versus 
storage time, C(t) = w(t)/w(0), is given by C(r) 2 = [s(t) 4 + p 2 ]/[s(t) 2 (l + p 2 )], where 
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FIG. 3: (a) Linear growth of w 2 with respect to the time duration of diffusion. The line is 



W — Wn 



ADt, with D = 10.8 cm 2 /s. Both D and wq were fit parameters (wq = 0.4 ~ 0.55 mm, 



varying for the different modes), (b) Self-similarity: cross-sections (cartesian or radially- weighted) 
at different times are congruent when plotted versus the scaled coordinate [y(r) = s{t)~ 1 I 2 v for 
the cartesian HG mode, f(r) = s 

( r )-l/2 r 

for the cylindrical LG modes]. The solid lines are the 
analytic forms of the elegant HG and LG modes. 



p = z/zr specifies the initial distance from the waist. As the waist radius increases during 
the diffusion, zr(t) increases and eventually crosses the observed plane (which z coordinate 
is constant). At this time, the maximal contraction C^ in = 2p/(l + p 2 ) is obtained, and 
thenceforth the beam expands indefinitely. 

Elegant modes, as opposed to the standard modes, are generally not self-similar under 
diffraction, and the shape of the beam depends on z/zr. Hence, even when z is held 
constant, the increase of zr during diffusion changes the transverse shape and breaks its self- 
similarity. However, the aforementioned 'common' modes, which are simultaneously elegant 
and standard, are self-similar under diffraction and thus also self-similar under diffusion 
even for z ^ 0. The HG modes of orders and 1 9[ and all LG modes with p = (the 
vortex modes) are such common modes. Far from the waist, at z > zr, these modes contract 
self-similarly. Fig. [5] presents the experimental result for a diffusing mode focused at a 



7 



1 



.75 

m 
S 

O c 
Q. 

o 
I- 

.25 



" 50 100 150 

Storage time [ps] 

FIG. 4: Decay of the total power in the retrieved images P(t) = f dxdy\ip{r)\ 2 . In the inset, P* 
compensates for the algebraic decay [s(r) known from Fig. [3] and collapses onto a single straight 
line in the semilog scale, yielding the homogenous decay rate 27 = (68 /^s) _1 . In the main graph, 
the dashed line is e _2rr and the solid lines are e~ 2 " ,T s(t)~ 2 ( n+1 \ demonstrating the faster decay 
of the higher-order modes. The difference between the LGo,i and the HGo,i, both having N = 1, 
is due to slightly different initial waist radii (wq). 




distance 8z R before the cell (p — 8), yielding C^ in ~ 1/4. 

Finally, we point out an intriguing instability phenomenon noticeable in Fig. [2(a) for the 
mode. During diffusion, the m = 2 vortex breaks-down into two vortices, probably of 
m=l. A decay of high-order vortices into lower-order ones has been seen also in optics, 
quantum fluids, and Bose-Einstein condensates. Here, several candidate mechanisms may 
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29], which evidently conserves the cross-section of 



be responsible for the imperfection 
the original vortex [Fig. Efb)]. 

In conclusion, we have shown that when elegant Gaussian modes are put through co- 
herent diffusion, their waist radius effectively expands. The total power in the field decays 
algebraically during the diffusion, even for the lowest-order mode, due to field interference 
effects. The complete set of elegant modes is self-similar at the waist, while far from the 
waist, self-similarity is found only for the subset of common modes. For the latter, a self- 
similar contraction may occur. The cross-sections at the waist of the standard and the 
elegant Gaussian modes form a complete set of self-similar modes for diffraction and diffu- 
sion, respectively. A similar set for the simultaneous process of diffraction and diffusion, as 
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FIG. 5: The self-similar contraction and subsequent expansion of the LGo,i mode upon diffusion. 
Contrary to the former experiments, here we refer to the original beam as the initial condition (left, 
taken off-resonantly) , since much contraction occurs already during the slow- light propagation in 
cell. The original mode has a radius of 700 ^m and a phase curvature of (250 mm) -1 . The 
minimal radius, of 340 /um, is obtained after slowing for 5 [is plus storage for 7 /is. The black line is 
calculated from C(t) 2 with no fit parameters. In these conditions, the power decays substantially 
faster than in the expansion experiments, and substantial noise is already apparent after 15 /xs of 
storage. 



occurring in the dynamics of slow light 21], is yet to be explored. 
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